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We study the evolution of maximally symmetric p-branes with a Sp-i ® IR* topology in flat 
expanding or collapsing homogeneous and isotropic universes with A*' + 1 dimensions (with > 3, 
p<N,Q<i<p). We find the corresponding equations of motion and compute new analytical 
solutions for the trajectories in phase space. For a constant Hubble parameter, H , and i = we 
show that all initially static solutions with a physical radius below a certain critical value, r°, are 
periodic while those with a larger initial radius become frozen in comoving coordinates at late times. 
We find a stationary solution with constant velocity and physical radius, rc, and compute the root 
mean square velocity of the periodic p-brane solutions and the corresponding (average) equation 
of state of the p-brane gas. We also investigate the p-brane dynamics for H 7^ constant in models 
where the evolution of the universe is driven by a perfect fiuid with constant equation of state 
parameter, w — Vp/pp, and show that a critical radius, rc, can still be defined for —1 < w < Wc 
with Wc ~ {2 — N)/N. We further show that for ui ~ TOc the critical radius is given approximately 
by rcH oc {wc — w)'^" with 7c — —1/2 {rcH 00 when w Wc). Finally, we discuss the impact 
that the large scale dynamics of the universe can have on the macroscopic evolution of very small 
loops. 
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I. INTRODUCTION 



The formation of topological defect networks is ex- 
pected at cosmological phase transitions in the early uni- 
verse P, |3l ■ The cosmological consequences of such phase 
transitions can be very diverse depending both on the 
type of defects formed and on the evolution of the uni- 
verse after the phase transition. For example, it is well 
known that the average domain wall density in an ex- 
panding universe grows with respect to the background 
density in the radiation and matter eras and consequently 
domain walls must necessarily be very light in order not 
to cornpletely dominate the energy density of the uni- 
verse [3|. On the other hand, this late time dominance 
is a property naturally associated with the dark energy 
and domain wall networks have been proposed as inter- 
esting dark energy candidates [3]. However, a domain 
wall network would need to have a very small character- 
istic length and velocity in order to be the dark energy. 
This requirement has recently been shown not to be ful- 
filled both in the simplest domain wall models and in 
more complex scenarios with junctions [HQ. 

There are other topological defects, such as cosmic 
strings, which lead to much less dramatic consequences. 
A cosmic string network, in the scaling regime, has an 
average density that is roughly proportional to the back- 
ground density. Consequently standard cosmic strings 
do not tend to dominate the energy density of the uni- 
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verse and they naturally generate a scale-invariant spec- 
trum of density perturbations on cosmological scales. Al- 
though recent cosmological data, in particular the cos- 
mic microwave background anisotropics, rules out cos- 
mic strings as the main source of density perturbations 
on large cosmological scales, a cosmic stringcontribution 
at a subdominant level is all but excluded [3, [sl . In fact, 
the interest in cosmic strings has recently been revived 
by fundamental work suggesting string production at a 
significant (but subdominant) level at the end of brane 
inflation 

Cosmic strings and other defects can also be formed 
during an inflationary era or in between several periods 
of inflation. It is thus crucial to understand their evo- 
lution in such regimes in order to quantify their ability 
to survive any inflationary period which might occur af- 
ter they form [l^, |Tl[ . The evolution of circular cosmic 
string loops and spherical domain walls in a flat FRW 
universe has previously been studied in [l3| (see also 
p^ ) where the existence of periodic solutions in a de 
Sitter universe has been demonstrated. Here, we gener- 
alize these results in particular by explicitly computing 
the phase space trajectories and determining the critical 
radius for spherical p-branes in space-times with an arbi- 
trary number of dimensions. We also study in detail the 
more general evolution of maximally symmetric p-branes 
with a Sp-i (E) IR* topology in expanding and collapsing 
flat Friedmann-Robertson- Walker (FRW) universes. 

The outline of this paper is as follows. In section II we 
derive the equations of motion for spherically symmetric 
p-branes in FRW universes with an arbitrary number of 
spatial dimensions. We generalize these results in Section 
HI to account for maximally symmetric p-branes with a 
Sp-i ® topology. In section IV we study in detail the 
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dynamics of p-branes in expanding and collapsing uni- 
verses with a constant Hubble parameter, H, computing, 
in particular, the phase space trajectories and the crit- 
ical radius, Tc, for spherically symmetric p-branes and 
discussing p-brane evolution in various limits. In section 
V we relax the constant H condition and consider the 
possibility of a universe dominated by a dense p-brane 
gas. We determine sufficient conditions for the cosmol- 
ogy to have an impact on the macroscopic dynamics of 
small spherically symmetric p-branes in section VI and 
we conclude in section VII. 



II. SPHERICAL SYMMETRIC p-BRANES 

The space-time trajectory of a p-brane in a 
dimensional space-time is described by the Nambu-Goto 
action 



(1) 



where Tp is the p-brane tension, ^ab = g^lvX^a^'^b, g^v are 
the metric components and — {(,"") represents the 
space-time trajectory of the p-brane (a, 6 = 0, . . . 

In a A^-|-l-dimensional Minkowski spacetime, the tra- 
jectory of a p-brane with spherical symmetry may be 
written, in hyperspherical coordinates, as 

x{t, (jJi,. . . , (pp-i) = q{t) [cos , sin cos 02, 
sin (pi sin (/)2 cos (/)3, . . . , 
sin(/)i . . . sin0p_2 cos0p_i, 
sin 01 . . .sin(?!)p_i,0, . . . , 0] , (2) 

where 0i, . . . , </)p_2 G [0,7r[ and 0p-i G [0, 27r[ and for 
simplicity the coordinate system was chosen in order to 
make the defect aligned with the first p-dimensions. The 
area of a p-dimensional spherically symmetric p-brane is 
the area of a p-dimensional hypersphere. 



where |g| is the physical radius and 

-r 072 + 1) ■ 



(3) 



(4) 



The energy of the p-brane is, then, proportional to S'p7, 
where v — dq/dt and 7 = (1 — v^)^^^^. It then follows 
from energy conservation that 



dR^ dv^ /P7^ 
dt ' dt ^ \ R 



(5) 



where R = \q\"/^^^ is the invariant radius and s = 
sign{-q). 

In a flat A^+l-dimcnsional Friedmann-Robcrtson- 
Walker Universe the energy of the p-brane is no longer 
conserved and, consequently, the value of R is expected 



to depend on time. Energy-momentum conservation for 
a planar p-brane implies that the momentum per comov- 
ing area (that is, the p-dimensional generalization in a 
A^+l-dimensional spacetime of the concept of the usual 
2-dimensional area in a 3-|-l-dimensional spacetime) is 
proportional to a^^ and consequently the velocity of a 
planar p-brane satisfies vj cx a"*^^^^^ which implies 



dv 
'dt 



(6) 



The equation of motion for a spherically symmetric p- 
brane has an additional curvature term (see Eqn. ([5])) 
and is then given by 



dv 



PI 



1/p. 



so that 



f=HR 
dt 



R 



1 



{p+l)vH 



P 



(7) 



(8) 



where R is now defined as R — ^^^P\r\ with r = aq. 
Note that, in this case, if iJ 7^ then R is no longer 
time- independent . 



III. 



) R p-BRANES 



Let us generalize the results of the previous section for 
maximally symmetric p-branes with a S'p_i(8)R' topology, 
with p < N and < i < p, by allowing a number i 
of dimensions to have no curvature and leaving p — i 
dimensions with spherical symmetry. In this case the 
area of the defect (per unit of i-dimensional area of the 
non-curved dimensions) is 



s;^ip^t + i)Cp-,+n\q\^ 



(9) 



and energy-momentum conservation leads to the follow- 
ing equations of motion 



dt ^ 



(p-i)7i/(P-*)s 



R 



(p + l)vH 



dR 
'dt 



HR 



p+ 1 



,(10) 
(11) 



where the so-called invariant radius, R, is now given by 
R = |g|7i/(p-i)a. Note that 



dr 
'dt 



Hr, 



V ^ ±\ 1- 



r \ 2(p-i) 

'rj 



(12) 
(13) 



These equations clearly show that there is an asymmetry 
between the collapse and the expansion of the p-brane. 
This asymmetry is clearly seen in Fig. 1 which is dis- 
cussed in detail in the following section. Note that If 
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the universe is expanding then there is a damping term 
which always contributes to decrease \v\. On the other 
hand, the curvature term may act to increase or decrease 
\v\ depending on whether the p-brane is coUapsing or ex- 
panding. It is also very interesting to realize that we may 
write the Hubble parameter as a function of the two loop 
parameters, R and v, as 

H^'-^(l-P-±lA~\ (14) 
dt \ p-i J 

What this equation clearly shows is that we could, in 
principle, infer the dynamics of the universe by looking 
at the dynamics of a single p-brane. 

Note that the equations of motion pil)l and (fTT|) are 
invariant with respect to the transformation g — > — g and 
t —f —t. This implies that the phase space trajectories in 
expanding and collapsing universes related by the trans- 
formation H ^ —H are identical. However, the trans- 
formation t —t implies that the direction in which the 
trajectories are travelled is reversed. 

IV. P-BRANE DYNAMICS (H=CONSTANT) 

We have determined the (r, v) trajectories in phase 
space in the case of spherically symmetric p-branes (i — 
0). Integrating Eqns. pil)) and ITTI) we have 

-yrPil + vHr) =k, (15) 

where 7 = (1 — w^)^^/^ and fc is a constant. A stationary 
solution of Eqns. pO)) and (fTTj) with fixed critical velocity 
{v = 0) and physical radius (r = 0), respectively Vc and 
Tc, is characterized by 

vl = H'rl = . (16) 

p + 1 

This solution describes a p-brane standing still against 
the Hubble expansion or contraction. If i = then 
= p/{p + 1) and we may easily find that value of 
k corresponding to the stationary solution is given by 
k = kc ^ |i/|-PpP/2(p -f l)-(P+i)/2. Trajectories with 
k > kc will approach the line defined hy v = —{Hr)^^ 
when |r| 00 (so that z; ^ in this limit). On the other 
hand all periodic trajectories have k < kc- Note that the 
reverse is not true since there is an infinite number of 
non-periodic trajectories with k < kc- This is clearly 
shown in Fig 1 where spherical domain walls trajectories 
with p = 2 and i = are represented in phase space. The 
trajectories with k < kc are represented by a solid line 
and the dashed lines represent trajectories with k > kc 
which start at the one of the critical points defined by 
V = —Hr = ±1 and end at (r = 00, z; = 0) in an expand- 
ing universe (or vice- versa in a collapsing universe) . Note 
that even the periodic trajectories, for k < kc, are asym- 
metric. This is a direct result of the asymmetry of the 
equations of motion (fT^ and (fT51) discussed in Section 
III. 
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FIG. 1: Representation of the trajectories of spherical do- 
main walls (p = 2, j = 0) in phase space. The critical points 
marked with a red dot represent the critical stationary so- 
lution defined in Eqn. [16] and the yellow dots correspond 
to the critical points defined by « = —Hr = ±1. The solid 
lines correspond to trajectories of spherical domain walls with 
k < kc (those with \r\ < describe periodic trajectories in 
phase space). The dashed lines correspond to trajectories 
with k > kc which start at the one of the critical points de- 
fined hj V = —Hr — ±1 and end at (r = 00, v = 0) in an 
expanding universe (or vice- versa in a collapsing universe). 

We may also compute a critical (initial) radius, r[!, cor- 
responding to an initially static p-brane solution (u° = 0) 
with k = kc (in an expanding universe it will asymptote 
to the stationary solution when t 00). This critical 
radius, is given by 

r^"^|i/|-V/2(p+l)-(P+i)/(2p). (17) 

In the case of circular cosmic string loops {p — 1) 
and spherical domain walls (p — 2) r° = \H\~^/2 and 
rj? — 2^/^3^^/'^\H\^^ respectively. We recover the critical 
radius, r^, for a circular cosmic string loop but the critical 
radius, r^, for a spherical domain wall is slightly smaller 
(by about 10%) than the approximate solution given in 
ref. [13] ■ For larger initially static spherically symmet- 
ric p-branes with \r'^\ > r|? the motion is not periodic 
and, if H > 0, the brane eventually freezes in comoving 
coordinates. Specifically, the p-brane will asymptotically 
behave as 

q = const. <=^ R cc a , v cc a^^ . (18) 

If -ff < then all solutions with |r"| > will asymptote 
the critical point defined hy v — —Hr — ±1 when t 00. 

Periodic solutions with |r°| < r° satisfy {d{lnR)/dt) = 
and consequently it follows from Eqn. (fTO|) with i — 
that 



4 




0.03 0.04 
Ht 



0.07 



FIG. 2: Time evolution of the invariant radius, R, of a domain 
wall with cylindrical symmetry with R'^H — 0.002 (solid line) 
and the predicted evolution of {R} (dashed line). 



where the brackets denote a time average over one period. 
Note that the r.m.s. velocity approaches unity in the 
p — > cxD limit. On the other hand, the energy density, 
Pp, in sub-critical isolated spherically symmetric p-branes 
should evolve as matter {pp oc a^^). Since the p-brane 
equation of state is w = Vnjpp — {{p + — p)/N 

then {v^) = p/{p + 1) [3, [l^ is required in order that 
w = 0. 

If i 7^ there are no longer periodic solutions. This 
can easily be seen in the limit where R\H\ 1. In 
this limit the Hubble damping term has a very small 
impact on the dynamics of the p-brane on timescales ^ R 
and consequently the dynamics of the p-brane is quasi 
periodic. Hence, from Eqn. we see that the evolution 
of the velocity is essentially the same for all branes with 
the same value of p — i. In particular 



p ~ I 
p — i + 1 



(20) 



which only depends on p ~ i. We may then write Eqn. 
(HU as 



dt 



1 



p — i + l 
p — i 



2 ' 2 
V V 

p — i 



(21) 



Taking the average of the r.h.s of (jlOp over one quasi 
period one obtains 



d{R) 
dt 



= H{R) 
= -H{R) 



1 - 



p - i + 1 ^ 
p — i 



p — I 



p - i + l 



so that 



(R) oc exp (oiiHt) 



(22) 
(23) 

(24) 



with ai = —ijij) — I + 1) . Hence, if i 7^ and the uni- 
verse is expanding then the p-brane radius shrinks over 
cosmological timescales (the opposite occurs if _ff < 0). 
In Fig. 2 we confirm the above results by comparing the 
evolution of the invariant radius, i?, computed numeri- 
cally with the analytical macroscopic solution given by 
Eqn. (|24p and we find a good agreement between them. 



V. P-BRANE DYNAMICS (H / CONSTANT) 

We shall now discuss the case in which the universe is 
expanding and the Hubble parameter is no longer time 
independent. For simplicity we shall assume that the 
dynamics of the universe is driven by a fluid with w ~ 
constant ^ -1 so that a cx \t\^ with /3 2/{N{w + 1)). 
We identify t — either with the big-bang (for w > — 1) 
or with the big rip (for w < —1). Ifw>— 1 the Hubble 
radius, H^^ increases with time but for —\<w<Wc 
with Wc — {2 — N)/N (so that (3 > 1) the comoving Hub- 
ble radius {H~^ / a) decreases with time and consequently 
it is still possible to find a critical radius, associated 
with a solution that has v Vc ^ Q when t ^ 00. By 
requiring that v = Q and v — Q asymptotically at late 
times and using Eqn. pU]) we find that 
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P — 1 p — i 



We see that in the /3 



(3 p+1- 
1 limit Vr ~* and 



rcH (X {(3 — oc {wc — wY' 



00 . 



(25) 
(26) 

(27) 



The value of the critical exponent, 7c is equal to —1/2 
and is independent of i, p and N. 

The critical radius defined above no longer exists for 
w < —1 or w > Wc- If w > Wc then (3 < 1 and con- 
sequently the comoving Hubble radius, H^^ /a, is an in- 
creasing function of time. This means that all p-branes 
(even those that are initially very large) will eventually 
come inside the Hubble sphere so that at late times, when 
r <C H^^, the p-brane will oscillate quasi-periodically. If 
w < ~l then P > 1 and the Hubble radius decreases 
with increasing physical time. Consequently, any initial 
quasi-periodically p-brane trajectory (with r <^ H~^) 
will eventually eventually freeze in comoving coordinates 
(with r > H^^). 

The analysis of the collapsing case is trivial an follows 
directly from the q — > ~q and t —t duality described 
in Section III. In this case, ii w < Wc then the comoving 
Hubble radius, |i?|^^/a, increases with time (if w < — 1 
then is also an increasing function of time). Hence, 

all p-branes will eventually come inside the horizon and 
will oscillate quasi-periodically when R\H\ becomes much 
smaller than unity. On the other hand, if w > Wc then the 
comoving Hubble radius, \H\~^/a, decreases with time 
and so all p-branes will eventually have a physical radius 
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much larger than and asymptotically the physical 

radius becomes \r\ (x a so that oc a~^P~^^\ Conse- 
quently, as the universe collapses (a — s- 0) the p-branes 
become ultra-relativistic while staying effectively frozen 
in comoving coordinates. This result has been demon- 
strated in ref. [IB, [13] for defect networks in 3+1 dimen- 
sions but remains valid when we consider the dynamics 
of defect networks in contracting FRW universes with an 
arbitrary number of spatial dimensions as long asm > Wc- 
An interesting case is that of a flat expanding universe 
dominated by a dense brane gas with average density Pp 
and average pressure Vp. The Einstein Equations for a 
N+1 FRW universe are given by 



N{N- 1) 



Pp ■ 



(29) 



where Gn+i is the dimensional Newton constant. If 
we assume a constant equation of state for the brane gas 
given hy w = Vp/ Pp = [[p + — p)/N (see previous 

section) then the dynamics of the evolution of the scale 
factor with physical time is given by a cx i'^ with 



/3 



iV(l-fw) N -p+{p+l){v^) 



(30) 



In order to accelerate the universe one needs /3 > 1 or 
equivalently w < Wc = {2 — N)/N and consequently 



N + p 



p+1 



(31) 



Here, the brackets denote an average over the brane gas 
network. We see that accelerated expansion is possible 



only if = p -f 1 and in that case one would need {v 
N^^ for inflation to take place. 



< 



VI. CAN COSMOLOGY HAVE AN IMPACT ON 
THE DYNAMICS OF SMALL COSMIC STRING 
LOOPS ? 

In Section IV we have shown that Hi then the cos- 
mology has an impact on p-brane dynamics over cosmo- 
logical timescales, even ii R\H\ <C 1. However, we have 
shown that if _ff = constant then the macroscopic evolu- 
tion of sub-critical spherically symmetric j»-branes (i = 0) 
is not influenced by the cosmology ((i?) = constant). In 
this section we ask whether or not there are cosmologi- 
cal models with a Hubble parameter, H{t), in which the 
macroscopic dynamics of small spherically symmetric p- 
branes can be affected by the large-scale evolution of the 
universe. 

If we consider the evolution of p-branes with R\H\ <C 1 
in a flat FRW universe then Eqns. (fTU)) and (fTT|) im- 
ply that the impact of the cosmology on brane dynam- 
ics is very small on a time scale < R. Also since, for 
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FIG. 3: Time evolution of the invariant radius, R, of a circular 
cosmic string loop for H{t) = Hq + Hi{l — 2v'^{t)) (solid line) 
and the expected evolution of the mean value of R (dashed 
line). We assumed have taken Ho = 2Hi and R°Ho = 0.002. 



H = constant, the trajectory of a spherically symmetric 
p-brane is periodic, the effect of the cosmology averages 
to zero on each period of brane motion. We will now 
assume that H is no longer a constant and is given by 
H{t) = Hq + AH{t) where Hq > is a constant and 
|Ai/| < Hq. For simplicity, in this section, we consider 
only the case of a circular cosmic string loop in 3-1-1 di- 
mensions (TV — 3, p — 1, i — 0). The generalization of 
our analysis for N, p and i arbitrary is straightforward. 
If we take AH ^ Hi{l - 2v'^{t + 9)), where v is the mi- 
croscopic velocity of the loop and is a constant phase, 
then Eqn. pil) becomes 



i? (l - 2v\t)) (Ho + Hi{l- 2v\t + 0))) . (32) 



dR 
'dt 

The evolution of the loop during one quasi period is 
hardly affected by the cosmological expansion and con- 
sequently the Minkowski space solution given by 



r = rocosi 

t 



(33) 
(34) 



is still a very good approximation, in the RH <C 1 limit, 
on time scales <C H^^. Taking into account that 1 — 
2v'^{t) = cos(2i) and an averaging the r.h.s. of Eqn. ([32|) 
over one period one obtains the equation that describes 
the macroscopic evolution of circular loop 



d{R) _ Hj_ 
dt ~ 2 



cos(26')(i?) , 



so that 



(R) oc exp {a2Hit) , 



(35) 



(36) 
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with OL2 = cos(2^?)/2. Wc clearly sec that the evolution 
of the universe may have an impact on the macroscopic 
evolution of cosmic string loops over cosmological time- 
scales, even if they are very small. To illustrate this we 
solved numerically the equations of motion for a cosmic 
string loop. In Fig. 3, we plot the results for the time evo- 
lution of the invariant radius of a loop with initial radius 
HPHq = 0.002 and for a Hubble parameter parametrized 
by by Hq = 2Hi and 9 = 0. As expected this cosmology 
has indeed an impact on the evolution of the invariant ra- 
dius, R, making its mean value increase after each period 
by the predicted amount. Of course, this is only valid in 
special cases, and in general such an effect is expected to 
be very small. 

VII. CONCLUSIONS 

In this paper we studied the dynamics of special p- 
brane configurations in expanding and collapsing FRW 
universes with an arbitrary number of dimensions. We 
have shown that, in the case of spherically symmetric 
p-brancs and a constant Hubble parameter, there are es- 
sentially two types of trajectories, closed periodic trajec- 
tories and open trajectories. We have obtained new an- 
alytical solutions for the trajectories in phase-space and 
computed the corresponding critical points as well as the 
root mean square velocity of the periodic p-brane solu- 



tions. We have shown that if i 7^ the solutions are 
no longer periodic with the p-brane radius evolving over 
cosmological timescales, even if R\H\ 1. 

We have also investigated the case where H ^ constant 
and we have found that for ~1 < w < Wc with Wc = 
(2 — N)/N a critical radius may still be defined. In the 
w ^ Wc limit we found that VcH oc {wc — w)^" with 
7c = —1/2. We also considered the special case of a flat 
FRW universe dominated by a dense brane gas and de- 
termined the required conditions for accelerated expan- 
sion to take place. Finally, we discussed the impact that 
the large scale dynamics of the universe can have on the 
macroscopic evolution of very small loops showing that 
there are situations in which the evolution of the universe 
as a whole may affect the macroscopic p-brane dynamics 
over cosmological timescales. 
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